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Abstract 

We extend the concept of quintessence to flat nonminimally coupled scalar- 
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tensor theories of gravity. By means of Noether's symmetries for the cosmo- 
logical pointlike Lagrangian C, it is possible to exhibit exact solutions for a 
class of models depending on a free parameter s. This parameter comes out 
in the relationship existing between the coupling F(tp) and the potential V(ip) 
because of such a symmetry for C. When inverse power-law potentials are 
taken into account, a whole family of exact solutions parametrized by such 
an s is proposed as a class of tracker fields, and some considerations are made 
about them. 

PACS number(s): 98.80.Cq, 98.80.-k, 98.80.Hw, 04.20.Jb 
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I. INTRODUCTION 



Recently, astronomical observations have indicated a strong evidence of an accelerated 
universe |T|-^|. Together with measurements of the cosmic microwave background and the 
mass power spectrum (see ||, for example), they suggest that a large amount of the energy 
density of the universe should have a negative pressure. A way of describing the missing com- 
ponent of the energy needed to reach, for instance, the critical energy density is quintessence 
Essentially, this is a spatially inhomogeneous and slowly evolving fraction of the to- 
tal energy. We may consider it as given by a scalar field (p slowly rolling down its potential 
V(ip) and such that — 1 < w v < 0, being w v = p v / 'p^, where p v and p v are, respectively, 
the pressure and the energy density of the scalar field. Actually, recent considerations f| fix 
the interval — 1 < w v < —0.6 as the more suitable one for such a scalar field to effectively 
represent quintessence. (A cosmological constant A, which mimics vacuum energy density, 
also produces a negative pressure, but this is such that Pa/ Pa = —1.) 

Within the scenario created by quintessence, there is a twofold problem. One of its 
two aspects is the so-called fine-tuning problem, based on the question why p^ appears to 
be so small with respect to typical particle physics scales. The other aspect, called cosmic 
coincidence ||, requires that the initial conditions have to be set precisely in order to explain 
why p^ and the matter energy density p m should appear of the same order today. This poses 
problems on the theoretical choice for the energy fraction which seems to be missing. 

More recently, a special form of quintessence has been introduced to avoid the coinci- 
dence problem. It is called tracker field |7|JT0|] and works like an attractor solution to the 
equations of motion, even if it is not really a fixed point, since it is time dependent and 
Pip/ ' Pm changes as (p evolves, leaving later cosmology independent of the early conditions. Of 
the kinds of potentials proposed for quintessence, two have been more studied for tracker 
solutions, namely, V(ip) = M 4+a ip~~ a and V(ip) = M A [exp(Mp/ip) — 1], where M and a > 
are free parameters, and Mp is the Planck mass. The family of tracker solutions is then 
parametrized by M, whose value can be fixed by the measured value of fl m today. Such 
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specific forms of potentials have been chosen because of their importance in particle physics 
models fTT|-|T4|| . Anyway, to our knowledge, Ratra and Peebles [15|;[16| were the first ones to 



study the influence of an inverse power-law potential in cosmology with a scalar field. 

A general study of specific features of tracker solutions has been made |7|,|10j , also leading 
to the introduction of the important function T, in order to fulfill the so-called tracker 
equation. Essentially, it is shown that "tracking behavior with wq < Wb occurs for any 
potential in which T = V"V / {V') 2 > 1 and is nearly constant" for any possible initial value 
of the scalar field Q. (Here, prime denotes derivative with respect to Q, and B indicates 
background.) As a consequence, once a certain potential has been assigned, the existence of 
tracking solutions can be tested without solving the equations of motion. 

Usually, the scalar field representing quintessence has been considered as minimally cou- 
pled to gravity and only more recently nonminimal coupling has been introduced in such a 
context fT7HT9|. In this connection, we consider interesting to refer again to scalar-tensor 
theories of gravity, in which the scalar field <p is nonminimally coupled to gravity and also 
inverse power-law potentials for ip have been studied, leading to exact solutions for ip(t) and 
for the scale factor a(t) of the universe (see |2DJ and references therein for a review on this 



topic). Our main purpose is to outline how a family of exact tracker solutions can be derived 
from a flat nonminimally coupled theory with inverse power-law potential, expliciting also 
how the tracker equation is fulfilled in such a way. 

In what follows, section 2 is devoted to a short review of some basic notions of flat 
nonminimally coupled theories, and section 3 selects a special class of solutions. In section 4 
we identify such a class as a family of tracker solutions, and in section 5 we draw conclusions. 

II. NONMINIMALLY COUPLED THEORIES 

As it is well known, one of the main reasons why nonminimally coupled scalar-tensor 
theories of gravity have received a special degree of attention is that they seem to play an 
important role in inflationary cosmology (see [2D], for instance). In the following, however, 
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we will limit ourselves to concentrate on what can be useful here, trying to draw a sort 
of narrow and straightforward path to our specific goal, namely deriving tracker solutions. 
First of all, we deal with flat (i.e., with the curvature scalar k — 0) models described by the 
action 

A = J d A x^ [F{<p)R + f<T^¥V - V(<p) + C m ] , (1) 

where g is the determinant of the metric g^, R the curvature scalar, semicolon indicates 
ccovariant derivative, and the functions F(ip) and V(ip) are not specified; C m is the La- 
grangian of an ordinary perfect fluid noncoupled to the scalar field (p. F(<p) expresses the 
nonminimal coupling of <p with gravity and is such that, when F = Fq = — 1/2 (using units 
such that 8irG — c — h — 1), action in Eq. (|l|) reduces to the usual one in the flat minimal 
coupling case. 

The field equations can be derived by varying Eq. (JTJ) with respect to and they can 
be written as 

fl _ rp _ _ rp(tot) (n\ 

vrfiv — -Lfiv— 2FUf) ^ V ' 

where G^ u = R^ u — \g^ v R is the Einstein tensor, and T^ u is a quantity related to the total 
stress-energy tensor 

rp{tot) _ rp{tp) _|_ J>(m) ^ 

Here, the tensor 



1 

T 

+2g^nF(ip)-2F{ip).^ (4) 



represents the scalar field source, while Tffl is the standard perfect fluid matter source (and 
□ is the usual d'Alembert operator). Varying with respect to ip, we get the Klein-Gordon 
equation ruling the dynamics of the scalar field p 

n ( p-RF'{<p) + V'(<p) = 0, (5) 



denoting the prime the derivative with respect to ip. It is possible to show that Eq. ([5]) 
is nothing else but the contracted Bianchi identity pl|,20|, which means that the effective 



stress-energy tensor introduced in Eq. (|2]) is a zero-divergence tensor, coherently with 
Einstein's theory of gravity [p^j . 

Fixing a homogeneous and isotropic (FRW) metric reduces the relevant variables to a 
and ip, i.e. the scale factor and the scalar field, each one a function of t only. As a matter 
of fact, field equations (H) can be reduced to two ordinary differential equations {k = 0) 



H* + !Lh+^- + ^ = 0, (6) 
F 6F 6F K ' 



H -AF~2\ H + 6Fj-2F + AF + 6F' (7) 

where dot indicates the time derivative, H = a/a, p m = w m p m is the equation of state for 
ordinary fluid matter, and p v = w v p v is the equation of state for the scalar field <p, having 
defined its pressure and energy density, respectively, as 

P V = W - V(cp) , (8) 



Pv = W + V(<p) • (9) 



This implies, thus, that 



^ - 2V(<p) 

Now, it is very interesting, and for us very important, to notice that Eq. (|5|) (rewritten 
in the FRW flat case) and Eq. (|7|) can be seen as the Euler-Lagrange equations of the point 
Lagrangian 

C = 6aa 2 F{cp) + 6dipa 2 F'(ip) + a 3 (p v + p m ) , (11) 

Eq. (^|) being equivalent to Ec = 0, where Ec is the energy. The configuration space is 
then given by Q = {a, ip} (the minisuperspace) and the tangent space by TQ = {a, a, ip, ip}, 
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being the coordinates associated with the pointlike Lagrangian £ just the scale factor a 
and the scalar field (p, with velocities a, (p. It is important to stress that from Eq. ( |TT|) 
we get that Ec = const.; making the homogeneous and isotropic limit of Einstein field 
equations (0) implies to choose such a constant equal to zero. According to Noether's 
theorem, the existence of a symmetry for the dynamics derived from £ involves a constant 
of motion. As a consequence, Noether symmetries in cosmology give the possibility to infer 
some transformations of variables which often lead to deduce exact cosmological solutions 
23| , |2^^0| . Such solutions, though obtained by means of a procedure suggested by the 



existence of this kind of symmetries, are actually independent of it and could also be got by 
chance, suitably choosing the right transformation of variables. That is, a Noether symmetry 
simply suggests that such a transformation should exist and gives an easy way to find it: once 
we have the right way to write down equations, we can easily solve them and get a solution. 
We have to stress, also, that, in order to verify the existence of a Noether symmetry, we 
find a way of assigning the functions F((p) and V(ip) for which a Noether symmetry exists, 
leading to remarkable results in many cases (see again |2D| and references therein for several 
examples). In the following, we will introduce a class of exact solutions which deserves a 
special attention in the context we are working in. 

III. A SPECIAL CLASS OF SOLUTIONS 

First of all, let us notice that we start from the pointlike Lagrangian (|TTJ) . It represents 
a whole class of theories, since a particular model is assigned by specifying F(ip) and V(<p). 
Examining the existence of Noether symmetries when the matter content is dust (p m = 0), 



as a first result P5| , |20| one deduces the relevant relation 

V(<p) = VoFfrjM , (12) 
where Vq > is an arbitrary constant assumed always positive, and 
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Therefore, the potential, through its exponent p(s), depends on the free parameter s. (The 
case s = —3/2 is degenerate and has been studied separately p6|,p0|. When s = —1, p(s) is 



zero and V(<p) becomes a constant; this situation has also been treated apart p0| .) 

The existence of such symmetries for Lagrangian ([0]) implies that we can find a differ- 
ential equation for F(tp) which has a general solution expressed by an elliptical integral of 
second kind, but we will limit our attention here to the particular solution 



F((p) = k ^ 



(14) 



where k < is an arbitrary constant. (Negative values of F(<p) are necessary to disregard 
repulsive gravity.) Then, Eq. (]i~4T) inserted into Eq. ([12]) gives 



V(<p) = [W (s) 



Ms) 



(15) 



It is possible to see that there is a whole family of exact solutions for the time evolutions 
of a and (p, and that they can be expressed as pO| 

(16) 



a(r) = Z(s)r r 



where 



V( T ) = C(s)t 



6s 



2s + 3 



X{s) 



m = a*) 



-2/[3p(a)] 



r = 



3/xW 



2s 2 + 9s + 6 



(17) 



(18) 



(19) 



s(s + 3) 

are parameters depending on s. Time r is, actually, a rescaled time. 

Apart from the cases showing pathologies in the solutions (s = and s = —3), which 
have to be discussed separately p5| |, it is important to notice that the right sign of the 
coupling, i.e. F(<p) < 0, implies 



2 < s < -1 . 
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(20) 



When s varies in such an interval, we have an infinite number of exact solutions of the forms 
given in Eqs. ([Li]) and (|iT|). Asymptotically, we have 

a(r) « r r , y?(r) « r - (2s+3)/s . (21) 

That is, depending on the values of s, the scale factor a(r) can have asymptotic Friedmann, 
power-law and pole-like behaviors. For instance, when \s\ 3> 0, it is a(r) ~ r 2 . But, in the 
range of values in Eq. (pCf), it can only be Friedmannian or power-law. 

As to the scalar field tp, it diverges for s < —3/2 and converges for s > —3/2. In what 
follows, the range of values anyway chosen for s will always be in the latter interval. 

IV. A FAMILY OF EXACT TRACKER SOLUTIONS 

Let us pose 

a = -2p(s)>0. (22) 
From the definition of p(s) in Eq. (pT3|), this implies 

-|<a<-l, (23) 

and Eq. (|l|) gives 

V(cp) = (v k - a / 2 )^ a , (24) 



with a always both positive and even. Substituting a into Eqs. (|Tq) , (^7|) , (Fq), and (|T9[) 
yields 

/2 "I- Q\ 4 /[ q ( 2 +«)] 2(2a 2 +9ct+6) 
(r) = ( — ^ ) T 3(a*+6a + 8) , (25) 

^(r) = (±±H) r 2 /^) . (26) 

On the other hand, from Eq. ( |ToD we see that, being 



*=w >0 (27) 

(where now dot indicates derivative with respect to rescaled time r), we can write 

x — 1 / \ 

u V = — . (28) 

Thus, it is clear that it is always w v > — 1. For a constant (p, i.e. a constant potential 
(which mimics a cosmological constant term), we should get w v — — 1. We have in general 

Cn 2 h^/ 2 

' s w =t kf >0 - (29) 

implying the constraint 

k a/2 > , (30) 

which is always true. 

One of the main requests that the scalar field has to satisfy, in order for it to be seen 
as a good tracker field, is that (p 2 < V(ip), i.e. that x < 1. Requiring that < x < 1 then 
poses the condition 

< < . (31) 

° 2(2 + a) y 1 

Of course, < x < 1 also implies that < 0, so that it is — 1 < w 9 < 0, namely what is 
needed for the scalar field to be interpreted as quintessence. Recently, it has been claimed 
that constraints from large-scale structure together with SNIa data imply w v < 0.6 with 
95% of confidence level ||, which forces the field of variation for x to be 

0<s< 0.25. (32) 

This, in turn, yields 

Constraints given by Eq. ([JT]) are also consistent with a direct calculation of from Eq. 
(|l0l). As a matter of fact, inserting into it the expression of Eq. (^) for y?(r), we get 
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That is, < implies Eq. fl3"ip. 

A further restriction can be found on the values of s, already such that Eq. ( |23|) holds, 
if one looks at equipartition at the end of inflation for inverse power-law potentials. This, 
being a even, constrains to a > 5 [|TUJ, so that 



3 21 

- < s < . 

2 16 



(35) 



In [10], in a minimal coupling regime, there were also introduced two important equations, 
the equation of motion and the tracker equation. As a matter of fact, the first one can be 
easily generalized to the nonminimal coupling situation, giving 

ldlnx 2F' 



V 

± v = 3 



1 



k(1 + w v 



1 + 



(36) 



6 d In a tpH 

where k = 8ttG/3 (we are changing now our units following the current litereture), Q v = 
Kp^/H 2 , and F' = 2k (f. (The ± signs, respectively, depend on whether V' > or V < 0.) 
F' — gives the minimal coupling case, and all cosmological solutions converge to the 
tracking solution, which is such that is nearly constant and less than wb (where B 
indicates background) , implying that 1 + = 0(1) and therefore (p 1 ~ Q^H 2 = np v , so 
that 



V 



H 





(37) 



This is referred to as the tracker equation [ |10|] . 

We can also introduce the function V = V"V/ (V) 2 , used in [|10| for a test on the tracker 
behavior. Let us notice, then, that we have at all times 



„ V"V 2p(s)-l , 1 
T = - — - = , r = 1 + - = const. 
(V) 2 2p(s) a 



(38) 



which is the major condition for a tracker behavior. It is T > 1 for s in the interval in Eq. 
([23]), being 
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4s + 3 



1 



2s + 3 



(39) 



6(s + l) * 6(s + l)' 
Of course, the range of values in Eq. ( j35|) for a > 5 is contained in the one in Eq. (p3|), still 
implying therefore T > 1. 

The tracker equation is then 

1 



1 + 



X 



2x 



0(1 + x) 2 n 



pi 



2 



k{1 + Wu 



X < 'x + 



6 + x + |f + if) 
(2/f 2 + if) 



d 

77^ 



12F' 
0ii 



(40) 



giving back what can be written for the minimal coupling case (when F' = 0). Here, 
x = din xj din a and x = (i 2 lnx/(ilna 2 . 

When an inverse power-law potential is considered, as shown in |27J for the minimal 
coupling case and in |Tj]] for a nonminimal coupling case (with F(ip) = £y? 2 /2, being £ a 
constant, but with a slightly different Lagrangian for ip), if the perfect fluid with ps oc 
a -3(i+w B ) dominates, so that a oc t 2 ^ 3 ^ 1+WB ^, then the following relation holds 

wbch — 2 



a + 2 



(41) 



On the other hand, if T is nearly constant, Eq. ( f4T|) in the minimal coupling regime implies 
fT0|| that a solution exists in which w v is also nearly constant and x, x, 'x become nearly 
zero. This also leads to Eq. (|4i~l). 

^From Eq. (|3^), relation in Eq. ([y]) can be written as 

2(r-i)(iw fl -i) 



w B - 



(42) 



i + 2(r-i) ' 

evidentiating the fact that T > 1 is equivalent to wb > w v in a matter dominated situation. 

In [JT7] it is found a solution for of the same type as in Eq. (|26|). Thus, the situation 
described therein is practically similar to ours, and we can import some of its considerations. 
Assuming that the universe is matter dominated, a oc 

t 2/[3/(i+w B )] implies tnat H = [2/3(1 + 
Wfljjt" 1 , and the Klein-Gordon equation @ can be written as 
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2 1 

(p + — y? 

1 +W B T 



1 + w B ) 



3(1 + w B ) 



- 1 



-aV k if 



a/2 - a -l 



0. 



(43) 



As demonstrated in ||17|| , this involves, for example, that we can take Eq. ([11]) as valid also 
in our context, even if we did not write Eq. ([H]) in such a way to clearly evidentiate that 
behavior when x,x, 'x are negligible and w 9 is constant. 

If we, only in a speculative way, compare Eq. (|34D , "exact" and always valid at any time, 
and Eq. (f|l|), "approximated" and valid only when p B ^> p v , we find 



4(2 + a 

This forces the constant k a ^ 2 to be 

ko a/2 = I 



a/2 a(l + w B ) T/ 4 + a(l - w B ) 



54 



(44) 



27V n 



2 2(2 + a) 



(45) 



and, according to Eq. (|5ID , leads to 

2 + a 
54 



< V < 



2 + a 
27 



(46) 



or, from Eq. 



to 



4(2 + a) 2 + a 

< Vo < 



135 



27 



(47) 



Also, substituting Eq. (0) directly into Eq. (|3lD and Eq. (p3|), respectively, gives 



or, equivalent ly, 



2 4 -3a 
w B < - , w B < — 

a ba 



2 4 
a < — , a < — 

w B bw B + 3 



(48) 



(49) 



Let us notice, then, that the second relation in Eq. fl4*8| ) always implies w B < (which is 
not good for ordinary matter), while the first relation in Eq. fl4"8] ) gives the right constraint 
w B < 1. Also, the second relation in Eq. (|49f) yields a < 4/3 for dust (w B = 0) and a < 1/2 
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for stiff matter (wb = 1), while the first relation in Eq. ( (49| ) lets us accept any value of a 
for dust, and gives a < 6 for radiation (wb = 1/3) or a < 2 for stiff matter. Finally, we can 
say that these considerations seem to imply that Eq. (^) limits too much the variability of 
x. 

Let us also notice that Eq. ([11]) could be read as 

2-(2 + a)K| , . 

w B ~ . (50) 

a 

Now, for ordinary matter (i.e., when < Wb < 1), it comes out (always disregarding the 
approximated equality) 

?—^<\wJ<^—, (51) 

2 + a ~ 1 ^' ~ 2 + a ' V ' 

implying (for a = 2) < w v < 1/2, or (for a = 4) < < 1/3, for instance. This 
means, then, that values w v < 0.6 are however possible, even if considerations made above 
let us understand that there may also be values of w v such that 0.6 < w v < 0. 

Of course, all these kinds of considerations have to be taken just as indicative, since Eq. 
(O) is not always valid and is usually read in one way: once we assign a specific value of wb, 



then it gives an approximated value of w^. For example, wb = gives ~ —2/ (2 + a) < 
for any a > 0, and from wb = 1/3 it is found w v ~ (a — 6)/[3(2 + a)} (giving the critical 
value a = 6, such that a > 6 implies positive and a > 6 implies negative w v ). 

On the other hand, anyway, taking the tracker condition in Eq. (|37D into account, we 
can immediately control its validity, being 

V a 1+w v % 



V if 
H 



T 



l+wz 



(52) 



r 1+ -s . (53) 
<P 

As a global consequence, we can assert that (^) is a good family of exact tracker solutions, 
parametrized by the constants Vq and /co, being 

M 4+Q = V ko a/2 . (54) 
14 



If we consider relations in Eq. fl44|) , apart from the contribution of the ordinary matter 
content, the parameter M depends only on s 

(2s+3)/[2(s+3)] 

(55) 



M 



As 



l l+w B ) + 6 



81(s + l)(2s + 3) 

The observational constraint £l v ~ 0.7 today, when ip rs O(Mp), implies that V(ip ~ Mp) ~ 
p mo , being p mo w 10~ 47 GeV 4 the current matter density This gives 

V(p) = M 4+ <V~ a « M 4+a Mp a w p mo , (56) 

so that M « (/j^Mp") 1 '^ > lGeV, in good comparison to particle physics scale, when 
a ^ 2, that is s ^ —1.2. (This shows again that, since such a value does not respect Eq. 
(|35|), constraining to a > 5 should be reconsidered.) 



V. CONCLUSIONS 

In the picture offered by observations in recent times, a good cosmological scenario 
needs an energy component with negative pressure. Of course, the simplest and most ex- 
treme candidate is a cosmological constant, but other softer proposals exist. Among these, 
quintessence has been advanced resuming, generalizing, and suitably readapting older ideas 
on cosmology with a scalar field. Initially thought of as a field acting like an attractor on 
other solutions, a more refined version of it has been proposed more recently, namely the 
tracker field. This kind of field tends to isotropize the universe at late times, nicely solving 
the coincidence problem. In that, it finds an ideal convergence with what is claimed by 
Wald's theorem on isotropization in homogeneous cosmologies with a positive cosmological 
constant. 

In this paper, we have shown that the requested features of such a field when V(ip) has 
an inverse power-law behavior can be obtained looking at a kind of exact solutions already 
present in the past literature on nonminimally coupled scalar-tensor theories of gravity (see 
|20| , for example). Even if it was deduced in a very peculiar way there, nevertheless it has 
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been reintroduced as, say, an ad hoc tracker solution only more recently. For these reasons, 
it has seemed interesting to connect the solution for ip to the solution for a, placing all the 
discussion in the context of a well developed theory. To be precise, then, a family of exact 
tracker solutions has been found here, depending on the values of a parameter s which is 
crucial for the relationship existing between the coupling F((p) and the potential V(ip) in 
our context. 

Such a relationship is a condition for a Noether symmetry to exist in the cosmological 
scenario involved by the model proposed through action in Eq. ([!]). This condition appears 
to be very important to us, because it seems to imply tracker fields in a very natural way, 



based on the apparent naturality of Noether symmetries in cosmology [^0 . 

Also, in the context of nonminimally coupled scalar-tensor theories of gravity, there is no 
dramatic difference between quintessence and cosmological constant proposals, introduced 
to solve the puzzle offered by recent observations. It can be shown that there is, in fact, an 
evident possibility to generalize Wald's theorem in order to get an asymptotic cosmo- 



logical constant in many nonminimal theories, without introducing it a priori (see [^] and 
references therein). As a matter of fact, for many choices of F(ip), a time dependent A-term 
can be defined, asymptotically approaching a constant (that is, a cosmological constant). 



(In |29j several examples are considered.) In this connection, it is noteworthy that some 



interesting comments have also been made |3(| on the fact that, with respect to the coin- 



cidence problem, an inverse power-law potential is not really different from a cosmological 
constant. 
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